FRW in cosmological self-creation theory: Hamiltonian approach 
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We use the Brans-Dicke theory from the framework of General Relativity (Einstein frame), but 

now the total energy momentum tensor fulfills the following condition ^ ^87rT'"^ M ' + T M "^'^ j = 

0. We take as a first model the flat FRW metric in the Hamilton- Jacobi scheme and we present 
the Lagrange-Charpit approach in order to find classical solutions. In the quantum scheme, once 
we determine the characteristic surfaces, the quantum solution is obtained. These two classes of 
solutions are found for all values of the barotropic parameter 7. 

PACS numbers: 04.60. Kz, 12.60.Jv, 98.80.Jk, 98.80.Qc 



I. INTRODUCTION 

The Scalar- Tensor theories have their origin in the 50's. Pascual Jordan was intrigued by the appearance of a new 
scalar field in Kaluza-Klein theories, especially in its possible role as a generalized gravitational constant. In this way 
appeared the Brans-Dicke theory, with the particularity that each one of the energy momentum tensor satisfy the 
covariant derivative 2], T^ 1 - 1 ;„ — 0, where i corresponds to the i-th ingredient of matter content. A few years latter 
(1982), appeared a new proposal by Barber known as self-creation cosmology (SCC) 0,0]. Since the original 
paper appeared in 1982, more and more authors have worked in the different versions of this theory in the 
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classical fashion. By instant, Singh and Sin gh 1C [ have studied Raychaudhary-type equations for perfect fluid in 
self-creation theory. Pimentel 11 1 and Soleng [12 



13] have studied in detail the cosmological solutions of Barbers self- 



ly ij 

20], £ 



14| . Venkatcswarhi and Reddy [3], Shri and Singh |l6l Il7l . Mohanty et al. [18} . Pradhan 



Sahu and Panigrahi [21], Venkateswarlu and Kumar [22[ are some of the authors who 
have studied various aspects of different cosmological models in self-creation theory, many of them usually adopt a 
particular ansatz for solve the Einstein field equations. These papers adapted the Brans Dicke theory to create mass 
out of the universes self contained scalar, gravitational and matter fields in simplest way. 

The gravitational theory must be a metric theory, because this is the easiest way to introduce the Equivalence 
Principle. However always it is possible to put additional terms to the metric tensor, the most obvious proposal is 
a scalar field <j>. Recently Chirde and Rahate [23j investigated spatially homogeneous isotropic Friedman-Robertson- 
Walker cosmological model with bulk viscosity and zero-mass scalar field in the framework of Barber's second self- 
creation theory and found classical solutions, it is the simplest way to work with this theory because only takes in 
account the energy-momentum tensor of usual matter and the scalar field <f>. 

This work is arranged as follow. In section II we present the method used in general way, where we take the Brans- 
Dicke Lagrangian density and we consider this in the self-creation theory. In section III, employing the flat FRW metric 
as a toy model with a barotropic perfect fluid, we found the energy density in dependence of the scalar factor and the 
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scalar field. In Section IV we construct the Lagrangian and Hamiltonian densities for the cosmological model under 
consideration and it is presented the classical solution using the Hamilton- Jacobi approach. The classical solutions 



are obtained using the Lagrange-Charpit method |2J-|26|, which structure is such that we find classical solutions for 
all values of the 7 parameter. In this same section we solve the corresponding Wheeler-DeWitt equation for this case. 
The section V we present the conclusions of the work. 

II. SELF CREATION COSMOLOGY IN GR 

The Lagrangian density in the Brans-Dicke theory is 

£[9, <t>] = V~9 (R<P - ^ ) 9' J ' U <i>, l l, < t > ,") + V^gLmatter, (1) 

where L matter = 16np, so making the corresponding variation to the scalar field and the tensor metric, the field 
equations in this theory become 

tf-^rt^ + ^n^-o, (2) 



R «0 _ lg^ R = ^ + ^ _ IgafyA^ + i ~ g^D 2 0) , (3) 

where the left side is the Einstein tensor, the first term in the right side is the corresponding energy momentum tensor 
of material mass coupled with the scalar field <f> . The second and third term corresponds at energy momentum 
tensor of the scalar field coupled also to <j> ■ Both equations are recombined given the relation 

U(f> = 4ttAT, (4) 

where A is a coupling constant to be determined from experiments. The measurements of the deflection of light 
restrict the value of coupling to |A| < 10 _1 . In the limit A — > 0, the Barber's second theory approaches the standard 
general relativity theory in every aspect. Ocf) — <j)\^ is the invariant D'Alembertian and T is the trace of the energy 
momentum tensor that describes all non gravitational and non scalar field matter and energy. Taking the trace of the 
equation §5§ and then substitute in (@}, we obtain the following wave equation to cj) 

comparing both equations Q and ([3]), we note that A = ( 3+ 2 2a ,) in the Brans-Dicke theory, where a; is a coupling 
constant. In the Barber's theory, the A parameter is a new coupling constant. 
Now the Einstein equation can be rewritten as 

QQ/3 _ ^rp Q/ 3(M) _j_ _J_rpQ^(0) _ rpQ^(T)^ /gN 



where 



In Brans-Dicke theory, each one energy momentum tensor satisfy the covariant derivative, T MI/ ; „ = 0. In the self 
creation theory, we introduce that the total energy momentum tensor is which satisfy the covariant derivative, 
T^ T ) .„ = 0, namely; 



T«^ T > . = 0, 



= 0, (7) 
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which imply that [Q a/3(T) ] = ^jr [Q Q/3(T) L where Q a ^ T ) = </>T a ^ T ). This equation is the master equation which 
gives the name of self creation theory, because the covariant derivative of this tensor have a source of the same tensor 
multiply by a function of the scalar field <f>. 

III. FRW IN SELF CREATION THEORY 

We apply the formalism using the geometry of Friedmann-Robertson- Walker 

dr 2 



ds" = -N 2 (t)dt^+A^(t) 



r 2 d6> 2 +r 2 sen 2 ((9)d(^ 



(8) 



1 — Kr 2 

where N is the lapse function, A is the scalar factor. Now we solve the equations ([3]), ((4]) and ([7]), with the aim to 
find solutions to density p(t), scalar factor A(t) and the scalar field <j>(t). 

Taking the transformation / = = 4- and using T Q/3 as a fluid perfect, first compute the classical field equations 
([3]), together with the barotropic equation of state P = jp, this equation become 

A'\ 2 A' <j>' lj (<i>'\ 2 p k 



the equation become 



AA 2 A6 2\ <h (b ' 6 A 2 y ' 



^ <t>" A .„ o \P rm 



The covariant equation (J7J or conservation equation to the total energy momentum take the following form 
A"4>' A'/0'\ 2 /0'\ 3 A'/0'\ 2 A' p S' p p> 

3 x0- 3w aUJ ~ a w - 24t i (1+7) r^ +8 7r 8l r 0, (12) 

thus, the system equations to be solved are (|9])-(fT2"]l. 

We write the term + ^r, in equation ([9]), using equations (fTQ)) and (JTTJ) , and after some algebra we have 

3 x " 3 + u (f ) 2 = 2?r [3A(37 ~ 1} ~ 2(1 + 37)1 i (13) 

Equation ([12")) can be rewritten as 

/A" A'0' /<A 2 W fU' P A' / P V 

|^ 3 X- 3 a7 + w UJ J? V a + Tj7"~ 247r(1 + 7) 0A - 8 "UJ ' (14) 

and using the equations (jTTJ) and (TTBl in (TT4")) . then we have the master equation for solve the energy density of the 
model, as 

2 [A(3 7 - 1)(3 + 2w) - 2(1 + 3 7 )] ^ - 24(1 + 7 )£^ - 8 = 0, (15) 



defining the function F = 4, we have 



oV Ln 



FA 3(i+ 7 ) i4[A(3 7 -i)(2w+3)-2(i+3 7 )] 



0. 
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i o = M 7 A- 3 ( 1+ ^V 3 , P = - ^ [2 - Aa ] , a = 3 + 2u;. (16) 



who solution is 



This equation is equivalent to General Relativity expression 4] with the addition the last factor representing the self 
creation cosmology. Note that for a photon gas 7 = |, so j8 = and equation (fTo) reduce to its General Relativity 
expression p = po (A/Aq ) ~ 4 which is consistent, because in the radiation epoch there was not interaction between 
photon and the scalar field. 

IV. LAGRANGIAN AND HAMILTONIAN DENSITIES IN SCC 

We will use the classical Hamilton- Jacobi approach to find solutions to (p, A, <fr) . The corresponding Lagrangian 
density using flSJ into (P) 

^ = ^ + ^-# + 1 ^N, (17) 



the momenta are (n g = 5$ 



N N N(f) 

A6A A 2 <j) A 2 A ujA 3 *]) 

A = 6(3 + 2 V 2 ( 3 ^ + - A "a), (18) 
* = 2(3+1^3 ( An A~2W, (19) 

when we write the canonical Lagrangian density ^canonical = n q q — NK, we obtain the corresponding Hamiltonian 
density as 

A -3 

H = — — -[-6(/) 2 Tll + wA 2 ni + 6A0n A n,i - 1927ra o A 6 0p], a = (3 + 2w). (20) 
12ao</> v 

A. Classical scheme: Hamilton-Jacobi equation 

Using the gauge N = 12a o 0A 3 , and the transformation in the momenta II q = where S is known as the 
superpotential function, with this, the Hamiltonian density is written as (we include the equation (fl6l) in this last 
equation) 

2 / no \ 2 



:2 



with Tj-y — 1927raoM 7 . In the following we obtain classical solutions for any value in the barotropic parameter p = 7/5 
using the Lagrange-Sharpit approach (^jj-Q. 

1. Lagrange-Sharpit approach 

For apply this method, we begin using the equation (|21[) . where we introduce the following transformations 

= e*, A = e n , (22) 



into (l2~Tj) . we have 



dS\ 2 /dS\ 2 /dS\ /dS 



6 'i Hm) +6 i i -'^ M ° e(S+1, *=°- (23) 



We start defining a function F($, p, q, S) = 0, where 

F($, n, p, q, S) = -6p 2 + wq 2 + 6pq - ^e^ 3 ^ 1 '^^* , (24) 

with the expressions P = §§ and g ~ ^ . We build the following set of equations 

dj = dO = dS = dp = dq = 

Fp dq pF p + qF q F$+pF s Fn + qFg ' 1 j 

calculating the derivatives we have 

F p = -12p + 6q, F q = 2wq+6p, F$ = ~(/3 + l)77 7 e~ 3(7 - 1)0 e (,3+1) *, 
Fn = 3( 7 - l^e- 3 ^- 1 )^ 1 '*, F s = 0, F t = 0, 

equation ([23)1 is read as 

d$ dn dS 



-12p + 6q 2wq + 6p p[-12p + 6q] + q[2wq + 6p] 

dq dp 



dt. (26) 



(27) 



-3( 7 - l)7 ?7 e- 3 ^- 1 ) f2 e^+ 1 )* (J3 + l)ri 1 e- 3 fr- 1 ) a eV>+ 1 ) i 

Choosing the equation 

dq dp 

-3(7- l^e- 3 ^- 1 )^^ 1 )* = (/3 + l)^ 7 e-3(7-i)n e (' 3 + 1 )* ' 

we find the following relation between the functions (p,q), 

{fi + l)q + 3(7 - l)p = c 7 = constant, (28) 

with this and from (|24p . substituting the term r] 7 e~ 3 ^^ 1 ' >il e^ +1 ' > '^ , we write an equation which involves to dependence 
of time, 

G8 + l)dp 



- dt, 
aip z + a 2 p + a 3 



where the solution is 



P(t) = ~ Va %~ 4aiS3 coth [c 7 V3a^t] , (29) 

a<x\ zai 

and the constants dj are 

ai ^ -6(/3 +1) 2 + 9^(1 - 7 ) 2 + 18(1 - 7 )(/3 + l), a 2 = 6c 7 [w(l-7)+^ + l], a 3 = wc 2 . 
From ([28)) we obtain the corresponding solution for the function q(t) as 

*> = Frr + TT>>- < 30 » 

Now we found the solution for the functions (p, q) , then we use the following equation in order to obtain the function 
*(*), 



being 



$(t) = a 4 t + a 5 Ln [Sinh(c 7 \/3a^t)] + $ ; 



where $o is an integration constant, and the constants 



£•4 



6c, 



3a 2 (3 7 + 2/3-l) 



/3 + 1 + 1 ai 
For obtain the solution for the function f2 we use the equation 



a 5 



6(3 7 + 2/3-l) 



ai 



dt 



2wq + 6p ' 



after integration, we have 



fi(t) = a6 1 — a,*? Ln [Sinh(cov / 3aot)] + f^o , 



where fio is an integration constant, and the constants 



1 



+ 1 



2w c 7 - — (l + /3 + w(l- 7 )) 
ai 



&7 



6(l + /3 + w(l-7)) 



(32) 



(33) 



(34) 



Now, taking in account the transformation law (|22|) . we have the solution for the scale factor A and the scalar field 
as 



A(t) = A e a6t Sinh- a7 (c 7 \/3aot), </>(t) = </> e a4t Sinh a5 (c 7 \/3a^t) . 
By introducing these solutions in the equation (fTTj) , the value of the constant c 7 is equal to 



(35) 



c 7 = v / 647raiM 7 A 2 </> 2 . (36) 
The corresponding value for all constants that appear in the classical calculation are presented in table [I] 
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Table I: All constants that appear in the classical solutions for various values ofj, with the following definitions £o 
li 



9+6u-A 2 a 2 ' 



36+24^ -A 2 a 2 . 



When we calculate the deceleration parameter 



AA 
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where the sign of q indicated whether the model inflates or not. The positive sign of q i.e. q > correspond to 
standard decelerating model, whereas the negative sign q < indicates acelerate expansion. 
Using (|35j) . this parameter have the following form 

q = —1 '■ j > (37) 

U/3ao^Cj cosh (c 7 ^/3aoi) — ag sinh (c 7 -\/3aoi)] 

and checking the table U over the possible value to the constants ao,a7 and c 7 , we observed that when uj > |, this 
deceleration parameter always is negative for any value in the barotropic parameter 7, so, the universe expand always 
in this theory. 



B. Quantum scheme 



Imposing the quantization condition and applying the Hamiltonian density (|20[) to the wave function we obtain 
the WDW equation for these models in the minisuperspace by the usual identification P q ^ = — id ql ± in ([20]) . with 
q' 1 = (A, <fr), and following Hartle and Hawking (27j we consider a semi-general factor ordering which gives (we include 
the equation (TTB1) in this last equation) (Also we could use the loop quantum cosmology in this approach [28|, [29[ ) 



H* = 



o 2 



d 



d 2 
dA 2 



<9A 



dA 



* = 0, 



(38) 



where Q and r are real constants that measures the ambiguity in the factor ordering between the scalar functions 
(A, <j)) and its corresponding momenta. 
Using the same transformation as the classical scheme, 



d 2 d 
6— 6(1 + 1) — 



d 2 d 
oj-tt^+uj(Q + 1) 



on 2 



on 



the equation 

d 2 



6 



dnd<P 



ry 7 e 



is read as 

3(7-x)n„08+i)« 



* = 0, 



(39) 



re-written this equation in the following form, 

d 2 *f d 2 ^ a 2 * 

6tt^t - 6- 



' d<f> 2 dfld<f> dn 2 
we can build the following a parameter 



9* 



■w(Q + 1)^ + ??7 e- 3 (^ 1 )°e^ +1 )** + 6(r + 1)^ . 



// 3ap 
6 : 



(40) 



which is used to construct a linear partial differential equation of first order of the characteristics surfaces for the 
functions ($, f2), given by 

89 + 86 

where is the envelope of the characteristic surfaces, then we have 

an de 



d$ 



' 



once we integrate this result gives us the following relation between these functions a ± $ — il = bo . 
Let f2 = <7 ± $ — bg, then introducing this in the equation (|39|) . we have 



8 2 V „, „,0* 
8<& 2 <9$ 



a 2 * _ ,,0* „ 9 2 * 

W ^+ w (Q + l)— -6- 



dn 2 v ^ 7 an sna* 

Assuming that the wavefunction is separable ^(0, $) = x/Ji(Q)ip2($>) we have 



e 3(7-l)bo e [0+l-3(7-l) CT ± ]*^ _ 0. 



6 d 2 V> 2 6(r + l)d-0 2 w d 2 7/>i w(Q + 1) d^i 



ip 2 d<& 2 



V-2 d$ Vi dO 2 



G 



1 di/'i 1 dip2 



dfl Vi d ^ ^2 d$ 



^ e 3(7-l)ho e [^+l-3(7-l)<r ± ]* = Q) 
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and considering the particular ansatz for solving this equation, ^^nr — Mo, which solution become xj}\ — bi e~^° n , 
and bi is a real constant. We choose only the minus sign, because we need to keep a decreasing wavefunction depending 
to scale factor. With this, the last equation is written as 



__ + b 2 _-(b 3 + b 4 e^ )^ = 0, 



(41) 



where the constants are defined as 



b 2 = /io - 1 - r, b 3 = [fx + 1 + Q] 



b 4 = ni e 3(i-i)*o 



b 5 =13 + 1- 3( 7 - l)a d 



equation that is similar to .''>( 



which solution is 



y = c 



y" + ay' + (be nx + c)y = 0, 



(42) 



where v — -^\/a 2 — 4:C, having the following relations between the parameters: 

x = $, a = b 2 , b= -6 4 , n = 65, c = -63, 
then substituting this relations and using the transformation equation <j> = e*, we have the following solution to ^2 



^2 = 



ciK 



v (k^'J + b 2 I„ (k^ 5 ) 



(43) 



where (K„, 1^) arc the modified Bcssel functions, /c = -^-j^—e <T 2 ' and £ = (1 + r — /-to) 2 + 2 " 3 A ' n [Mo + 1 + Q]- Then, 
the wavefunction have the following form 

$„ = A" M (j, 111 ^ [CiK„ + C 2 L • (44) 

The corresponding value for all constants that appear in the quantum calculation are presented in table [II] 
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2 + ,/' ,1/3 e- 60 

1+2<t± V 6 


2 /3T 

Aa n V 6 


2 /V-l -3b 
i>5 V 6 



Table II: j4W constants that appear in the quantum solutions for various values of 7. i/ere £ = (1 + r — /io) 2 + 2 " 3 M0 [/io + 1 + Q] 



V. CONCLUSIONS 



In this paper we have investigated the flat FRW cosmological model of the universe in the framework of Barber's second 
self-creation theory since of point to view of hamiltonian systems, where the classical solutions were found under the Hamilton- 
Jacobi approach combined with the Lagrange- Charpit method for all values of the barotropic parameter in a perfect fluid. By 
means of the deceleration parameter we obtain that the universe in this theory always suffers an expansion in any epoch of 
the universe. In the quantum behavior, the solutions were found using the curves characteristics method applied to partial 
differential equation of second degree, obtaining the solutions, as in the classical scheme, for all values in the 7 parameter. 



9 



Acknowledgments 

This work was supported in part by DAIP (2011-2012), Promep UGTO-CA-3 and CONACyT 167335, 179881 grants. 
JMR was supported by Promep grant ITESJOCO-001. Many calculations where done by Symbolic Program REDUCE 3.8. 
This work is part of the collaboration within the Advanced Institute of Cosmology, and Red PROMEP: Gravitation and 
Mathematical Physics under project Quantum aspects of gravity in cosmological models, phenomenology and geometry of 
space-time. 



[1] M. S. Berman and F. M. Gomide, Nuovo Cim. 74B, 182 (1983). 

[2] S. Weinberg, Gravitation and Cosmology: Principles and Aplications of the General Theory of Relativity, John Wiley and 

Sons, Inc. New York London Sydney Toronto, (1972). 
[3] G. A. Barber, Gen. Rel. And Grav. 14, 117 (1982). 

[4] G. A. Barber, A New Self-Creation Cosmology A semi-metric theory of gravitation Astrophysics and Space Science 282, 
683-730 (2002). 

[5] G. A. Barber, Self-Creation Cosmology ArXiv:1009.5862v2 (2010). 
[6] C. P. Singh, and S. Kumar, Int. J. Mod. Phys. D 15, 419 (2006). 

[7] C. P. Singh and, S. Kumar, Astrophys Space Sci Bianchi type-II space-times with constant deceleration parameter in self 
creation cosmology 310, 31 (2007). 

[8] C. P. Singh, S. Ram, and M. Zeyauddin, Astrophys. Space Sci. 315, 181 (2008). 

[9] M. K. Singh, M. K. Verma, and Shri Ram, Adv. Studies Theor. Phys. 6, 117-127 (2012). 
[10] T. Singh, Astrophys. Space Sci. 102, 67 (1984). 
[11] L. O. Pimentel, Astrophys. Space Sci. 116, 395 (1985). 
[12] H. H. Soleng, Astrophys. Space Sci. 138, 19 (1987a). 
[13] H. H. Soleng, Astrophys. Space Sci. 139, 13 (1987b). 
[14] D.R.K. Reddy, Astrophys. Space Sci. 133, 189 (1987). 
[15] R. Venkateswarlu, D.R.K. Reddy, Astrophys. Space Sci. 168, 193 (1990). 
[16] R. Shri, C. P. Singh, Astrophys. Space Sci. 257, 123 (1998a). 
[17] R. Shri, C. P. Singh, Astrophys. Space Sci. 257, 287 (1998b). 

[18] G. Mohanty, B. Mishra, and R. Das, Bull. Inst. Math. Acad. Sin. (ROC) 28, 43 (2000). 

[19] A. Pradhan, A. K. Vishwakarma, Int. J. Mod. Phys. D 11, 1195 (2002a). 

[20] A. Pradhan, A. K. Vishwakarma, Indian J. Pure Appl. Math. 33, 1239 (2002b). 

[21] R. C. Sahu, U. K. Panigrahi, Astrophys. Space Sci. 288, 601 (2003). 

[22] R. Venkateswarlu, P. K. Kumar, Astrophys. Space Sci. 301, 73 (2006). 

[23] V. R. Chirde, P. N. Rahate, Int J Theor Phys. 51, 2262-2271 (2012). 

[24] L. Elsgoltz, Ecuaciones Diferenciales y cdlculo variacional, Edit. Mir Moscu, (1969). 

[25] M. Delgado, SI AM Rev. 39, 298 (1997). The Lagrange- Charpit method. 

[26] G. Lopez, Partial Differential equations of first order and their applications to physics, World Scientific Pub. (1999). 

[27] J. Hartle and S. W. Hawking, Phys. Rev. D 28, 2960 (1983). 

[28] Xiangdong Zhanga, and Yongge Ma, Loop quantum Brans-Dicke cosmology, arXiv:1211.4183 . 

[29] Xiangdong Zhanga, and Yongge Ma, Loop quantum Brans-Dicke cosmology, J. Phys.: Conf. Ser. 360 012055 (2012). 
[30] Andrei D. Polyanin and Valentin F. Zaitesev, Exact Solutions for Ordinary Differential Equations, Chapman and Hall/CRC, 
(2003). 



